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The DGP brane world model allows us to get the observed late time acceleration via modified 
gravity, without the need for a "dark energy" field. This can then be generalised by the inclusion 
of high energy terms, in the form of a Gauss-Bonnet bulk. This is the basis of the Gauss-Bonnet- 
Induced- Gravity (GBIG) model explored here with both early and late time modifications to the 
cosmological evolution. Recently the simplest GBIG models (Minkowski bulk and no brane tension) 
have been analysed. Two of the three possible branches in these models start with a finite density 
"Big-Bang" and with late time acceleration. Here we present a comprehensive analysis of more 
general models where we include a bulk cosmological constant and brane tension. We show that by 
including these factors it is possible to have late time phantom behaviour. 



1. INTRODUCTION 

There are some intriguing problems with the standard 
ACDM cosmology that may hint at new physics in order 
to answer them. One of the problems is the requirement 
of adding an additional "dark energy" field in order to ex- 
plain the acceleration that the universe is now experienc- 
ing. An alternative way of explaining this phenomenon 
is to modify gravity at large scales. This should not be 
done in an ad hoc way but requires a consistent and co- 
variant physical grounding. The brane model put for- 
ward by Dvali, Gabadadze and Porrati (DGP) Q, and 
generalised to cosmological branes by Deffayet has 
such a grounding. 

Brane models are inspired by the discovery of D-branes 
within string theories. The key property is that matter 
is confined to the brane. The DGP model has a large 
scale/low energy effect of causing the expansion rate of 
the universe to accelerate. This is achieved via the addi- 
tion of an Induced- Gravity (IG) term to the gravitational 
action. The IG term is produced by the quantum inter- 
action between the matter confined on the brane and the 
bulk gravity. It has also been shown that the IG term 
can be obtained, in some string theory models, from the 
presence of Gauss-Bonnet gravity in the bulk Q . 

One of the problems with the DGP model is the fact 
that it is in some sense "unbalanced" . This is because it 
docs not modify the short-distancc/high-cnergy regime 
which we would expect from a string and quantum moti- 
vated theory. Brane models with Gauss-Bonnet (GB) Q 
bulk gravity on the other hand modify the high-energy 
regime like the original Randall and Sundrum (RS) mod- 
els did. The RS models have early times described by 
5D gravity because the bulk is warped. This warping 
localises gravity to the brane in the low energy regime 
while in the high energy era gravity leaks off the brane 
and "sees" the bulk and thus acts 5D. The DGP model 
has 5D behaviour once a certain length scale (the cross- 
over scale) has been reached. At this scale gravity starts 
to leak off the brane and into the bulk thus causing grav- 
ity to become strongly 5D. If we try warping the bulk in 
the DGP model with the aim of achieving both an early- 



time and a late-time modification, we ultimately fail. It 
has been shown 0| that the DGP model in a warped bulk 
either has 4D gravity behaviour at all scales or has 4D 
gravity at short and long scales with 5D behaviour at 
intermediate scales. 

We include GB gravity in the bulk with IG on the 
brane (GBIG) 0, If with the goal of obtaining a model 
that is "balanced", i.e. a model that gives us both UV 
and IR modifications. 

In Ref. we looked at the simplest generalization of 
the DGP model. We saw that this model provides us with 
both early and late time modification of GR. We found 
that by including both the GB and IG terms the "big- 
bang" singularity can have a finite density. In this work 
we systematically investigate the other solutions allowed 
in the model, i.e. we no longer restrict ourselves to a 
Minkowski bulk with zero brane tension. 



2. FIELD EQUATIONS 

The general gravitational action contains the Gauss- 
Bonnet (GB) term in the bulk and the Induced Gravity 
(IG) term on the (Z2 symmetric) brane: 
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where r > is the Induced Gravity 'cross-over' scale 
(r = k| / K4 consistent with Ref. 0, ) j is the Gauss- 
Bonnet coupling constant and A is the brane tension. 
A negative value of a would require the string coupling 
constant to be imaginary. We could still use the GB 
gravity in a classical sense but we would lose our string 
motivation. We have the standard conservation equation: 
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Using these variables we can write the two solutions for 
$ as: 
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The bulk cosmological constant gives us an upper 
bound on the GB coupling constant a.. Equation ((S} 
gives us: 
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For an RS (a — > 0) limit we take the minus branch. 
A5 > would imply £^ < and thus we must have A5 < 
and: 



FIG. 1; The effective bulk cosmological constant 
tion of X- 



• as a func- 



so the matter that is on the brane does not exchange en- 
ergy with the bulk and acts as a perfect fluid as in normal 
Friedmann solutions. The general Friedmann equation 
was obtained in Ref. 0- Here we assume that the bulk 
black hole mass is zero and the brane is spatially flat, 
but wc allow non-zero brane tension and bulk curvature. 
The general form of the Friedmann equation is then: 
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where $ is a solution to: 
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The bulk cosmological constant is given by, Ref. [li 
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Equations (0J and (0 give us two solutions for <&: 
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In dimcnsionless form this is given by: 
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Maintaining an RS limit would also rule out the 02 
branch. Wc would be restricted to the solutions lying 
along the line in the top left quadrant in Fig. We are 
interested in the whole range of the model so wc include 
the plus branch in Eq. ©. We assume A5 < therefore 
our constraint on a is given by: 
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If we take x — then we have no bound on 7 (apart from 
beingpositive and real). This is the case considered in 
Ref. H. 

In Fig. ^ we have the two (j) solutions plotted as func- 
tions of X- The two solutions with (/) = both live in a 
Minkowski bulk, all the rest live in an AdS bulk. Note 
that one of these AdS solutions {x = 0, (f> = —4/87) has 
A5 = but $ = — l/2a acting as an effective cosmolog- 
ical constant 0- We see that for any allowed value of 
(j> we can be on either of the two branches. This means 
that we need not consider the £md (j)2 solutions to the 
Friedmann equation separately. We therefore consider 
Eq. (|12|l in terms of (j). We define the maximum value 
of 7, for a particular value of (j>, that's allowed by the 
constraint equation: 



We work with the Friedmann equation in dimcnsionless 
form by defining the following variables (with r > 0): 
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The dimcnsionless Friedmann equation is: 



1 + 7 + 
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{h'-<P) = [h'-if, + a)]\ (14) 
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with the conservation equation now given by: where ' = d/ dr and h = a' /a. The Raychaudhuri and 

acceleration equations are given by: 

l^i' + 'ih{l + w)n^Q, (15) 



and: 
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FIG. 2: Solutions of the Friedmann equation {h vs /i) with 
negative brane tension {a = —2) in a Minkowski bulk (0 = 0) 
with 7 = 1/20. The curves are independent of the equation 
of state w. The arrows indicate the direction of proper time 
on the brane. 



3. FRIEDMANN EQUATION SOLUTIONS 



FIG. 3: Solutions of the Friedmann equation with positive 
brane tension (cr = 2) in a Minkowski bulk (0 = 0) with 
7 = 1/20. 



and positive brane tensions there are three solutions (this 
is not always true as we will show later) denoted GBIGl- 
3. There are four points of interest in Figs. [2 and O these 



The model we considered in Ref. is the Minkowski 
bulk limit (x = 0) of the 0i case, which we now see to 
be equivalent to x = ^ in the 02 case. The results we 
present below are in terms of and thus include both 
01 and 02 . We will consider the effect of including brane 
tension in a Minkowski bulk before we look at the AdS 
bulk cases. The Minkowski bulk case is given by — 0, 
since Eqs. (jSJ and © imply A5 = 0. 



3.1. Minkowski bulk (0 = 0) with brane tension 



• {fii, hi) and {fic, hi): The initial density for GBIGl- 
2 (fii) and the final density for GBIGl and 3 (no), 
found by considering dfj,/d{h^) = 0, are given by: 
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where the plus sign is for /i; and the negative sign 
is for /io. The Hubble rates for these two densities 
are given by: 



In Figs. 13 and |3| we can see the solutions to the Fried- 
mann equation in a Minkowski bulk. For both negative 
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where the sign convention is the same as above. 
The points (/Uj^ei ^i,e) have h' = |oo|. Therefore cos- 
mologies that evolve to /io, end in a "quiescent" 
(finite density) future singularity. This singularity 
is of type 2 in the notation of Ref. ■ 

(/Xi,0): This is the density at which GBIG3 "loi- 
ters". The point was found by considering h = Q 
and is given by : 



(20) 



In the case considered in Ref. we had cr = so 
III = 0. We can show that for a < GBIG3 will 
not collapse but will loiter at a density of /i = —a 
before evolving towards (/Xe, he), by considering h[. 
At ^ = -cr, /i = we get /ij = from Eq. ((TE|l . 
In a standard expanding or collapsing cosmology 
/i' < at all times. The evolution for a radiation 
dominated universe can be seen in Fig. ^ A dust 
dominated universe spends longer at fii . In Ref. 
they consider a loitering braneworld model. An im- 
portant difference between the two models is that 
in Ref. [l^ they require negative dark radiation, 
i.e. a naked singularity in the bulk or a dc Sitter 
bulk. Also in our model the Hubble rate at the loi- 
tering phase is exactly zero. The time spent at this 
point is only dependent on w and the density of 
the loitering phase is only dependent on the brane 
tension. 

(0, hoo)- This is the asymptotic value of the Hubble 
rate as ^ ^ 0. For the value of the parameters used 
in Fig. 13 GBIG2 is the only case with this limit. 
In general any of the models can end in a similar 
state (Fig. O . The different values of hoo are given 
by the solutions to the cubic: 
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In the case we considered in Ref. (o" = 0) the last 
term in the above cubic is zero. Therefore we have 
hoo = for GBIG3 while the solutions for GBIGl-2 
are given by: 
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where the minus sign corresponds to GBIGl and 
the plus sign to GBIG2. This is the only case where 
we can write simple analytic solutions as in all other 
cases we have to solve the cubic. 

The two Hubble rates, hi and he, are independent of 
the brane tension which simply shifts the ^ = axis. 
Therefore /io only corresponds to a physical (positive) 
energy density when a < < 0, see Eq. (|23|l . The effect 
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FIG. 4: Plots of /i vs r and h vs t for GBIG3 in a Minkowski 
bulk {(j) — 0) with negative brane tension {a — —1). We 
see that GBIG3 loiters around fi = ni — ~a, h = before 
evolving towards a vacuum de Sitter solution (/ic < 0). (Here 
7 = 1/20 and w = 1/3.) 



of the brane tension on the /i = axis is illustrated in 
Fig. El 

We showed in Ref. that for = = cr we require 7 < 
1/16 for GBIGl-2 solutions to exist within the positive 
energy density region (if 7 = 1/16 GBIGl-2 reduce to 
the same vacuum de Sitter universe). If we have some 
non-zero brane tension this constraint is modified. The 
maximum value of 7, for GBIGl-2 to exist with positive 
energy density, as a function of cr can be seen in Fig. 
We have defined new quantities, Co^i and cr/, for which 
Mc.i.i — 0. We can see from Eq. (|20|l that cr/ = 0. cto^i are 
given in terms of 7 by: 



1 - 187± (1 - 127)3/2 
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FIG. 5: Solutions of the Friedmann equation {h vs /i) in a 
Minkowski bulk (0 = 0) with 7 — 1/20. The vertical lines 
represent the /i = axis for the labeled brane tensions. 
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FIG. 6: The (it, 7) plane for solutions in a Minkowski (0 = 0) 
bulk. The short dotted horizontal hue is 7in = 1/12. GBIGl-2 
exist with positive energy density in regions III, IV and V. 



The maximum value of 7 (7ni) for GBIGl-2 to exist 
(i.e. for Eqs. H18() and H19() to have real solutions) is 
7,„ = 1/12 . Actually for 7 = 1/12, GBIG2 exists but 
GBIGl is lost. This is since /i; = /i^ at 7 = 1/12. The 
point hi = he is now a point of inflection and Eq. H19|) 
is still valid and therefore \h'\ — 00. For values of 7 > 
1/12, hi^c and /ii.o become complex, and it is no longer 
a point of inflection. Therefore \h'\ ^ 00 and GBIG3 
can continue its evolution through this point to h^o- In 
Fig.[7|we show results for h and h! for three values of 7, 
7 = 1/12 - 0.01, 1/12, 1/12 + 0.01. In the top plot we 
see that the solution denoted (1) has GBIGl-3 present 
as 7 < 1/12, we have two real and different values for 
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FIG. 7: The top plot \& h n& [i for three different solutions, 
all with = 0. The bottom plot shows h! vs for the same 
solutions. (1):7 = 1/12 - 0.01, a = -3. (2)17 = 1/12, a = 
-2.5. (3):7 = 1/12 + 0.01, a = -2. Different brane tensions 
are used purely for clarity. 



hi and h^,. This solution in the bottom plot has three 
parts (GBIGl has negative values not shown in Fig.[7J), 
GBIG2 is the dotted solution and the solid is GBIG3. As 
(1) approaches hi^^ \h!\ — * 00. Solution (2) has 7 = 1/12, 
in the top plot we see that GBIGl now no longer exists as 
hi = he. In the bottom plot GBIG3 (solid) and GBIG2 
(dotted) solutions both go to /i' = 00 at the point hi = h^. 
Solution (3) has 7 > 1/12. Now hi^c does not represent 
d^/dh = 0, therefore h' stays well behaved throughout 
the evolution (bottom plot). GBIG3 moves seamlessly 
onto GBIG2 making a single solution. As GBIG2 super- 
accelerates this new solution will show late time phantom 
like behaviour {h' > ^ w < —1). 

Solutions that lie along the ai line in Fig. |H| are those 
considered in . Regions I and II in Fig. |3 extend up 
to 7m = 00. Solutions in each region of Fig. Elhave the 
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FIG. 8: /loo for solutions in a Minkowski {(j> — 0) bulk. The 
thin-dark line has a = —2, thin-light line has a = —1.2, thick- 
dark lines have a = —1/2 and the thick- light lines have a = 2. 



following properties: 

• / : cro(7m) < cr < 0, 7 > 7i and a < CTc(7m), 7 > 
7m. GBIGl-2 do not exist. GBIG3 loiters at 
before evolving to a vacuum de Sitter universe. 

• // : (T > 0, 7 > 7i. GBIGl-2 do not exist. GBIG3 
evolves to a vacuum de Sitter universe. 

• /// : cr < (To, 7 < 7m- GBIGl will evolve to 
iHe^K)- GBIG2 evolves to {0,h^). GBIG3 loiters 
at m before evolving towards (/Xe,/ie)- When a = 
CTe GBIGl and 3 evolve to (0, /le). When 7 = 7m 
GBIGl ceases to exist {hi = he). 

• IV : cTe < cr < 0, 7 < 7i. GBIGl-2 both evolve 
to vacuum de Sitter states. GBIG3 loiters before 
evolving to a vacuum de Sitter universe. Each vac- 
uum de Sitter state has a different value of hca- 
When 7 = 7i GBIGl-2 live at (0, h). 

• V : cr > 0, 7 < 7i. GBIGl-3 all end in vacuum de 
Sitter states. With 7 = 7i GBIGl-2 live at (0, h). 
When cr = GBIG3 ends in a Minkowski state. 

In Fig.|Slwe can see how solutions in each of the regions 
mentioned above affect h^o. The thin-dark line (cr = —2) 
lies in III and I in Fig. Only GBIG2 exists for < 
7 < 1/12. For 7 > 1/12 we've seen that GBIG3 and 
GBIG2 connect to make a single solution, which is why 
the line is continuous through 7 = 1/12. The thin-light 
line (cr — —1.2) has an interesting feature due to the 
solutions lying on a line that cuts through both region 
III and IV as well as I. Both the thick lines have GBIGl- 
3 ending at /loo- The thick-dark solution loiters before 
evolving to hoc , but this has no effect upon the results in 
Fig.il 
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FIG. 9: Solutions of the Friedmann equation (h vs [i) with 
negative brane tension (cr = —2) in an AdS bulk (<^ = —0.1) 
with 7 = 1/20. The curves are independent of the equation 
of state w. The arrows indicate the direction of proper time 
on the brane. 



3.2. AdS bulk (<?!) 7^ 0) with brane tension 

When = 0, we have A5 7^ 0, with one exception: the 
02 solution with x = 0, i.e. 02 = —4/37, has A5 = 0, but 
the bulk is AdS, Ref. %. For x > 0, A5 7^ 0. Thus in all 
cases, 7^ implies an AdS bulk. 

When we allow the bulk to be warped (0 7^ 0) we open 
up another possible solution, denoted GBIG4. There is 
a maximum value of for which GBIG4 can exist as 
we shall see later. The nature of GBIG3 is also changed. 
These solutions, for a negative brane tension, can be seen 
in Fig. El Brane tension affects the solutions in the same 
manner as in = case. The qualitative effect of the 
warped bulk is to make GBIG3 collapse and to introduce 
the new bouncing branch GBIG4. This is due to /i = 
now giving two solutions: 



Mc,b = ±V^(2 + 70)-'7, (24) 

where the plus sign is for the GBIG3 collapse (/ic) and 
the minus for the GBIG4 bounce (/.tb). Effectively the 
loitering point in the Minkowski bulk is split into the 
max/min densities of the bouncing/collapsing cosmolo- 
gies of GBIG3-4. The other points in Fig.^arc modified 
by the warped bulk; they are now given by: 



1 - I87 + 277V ± (1 - 127 - 187^0)3/2 

^'■^^ si? 

(25) 

with the plus sign for [i\ and the negative sign for [if,. 
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The Hubble rates hi^c are given by: 
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/li o — \/l - 67 + 9-f^<t) ± VT~127 - 1872^, (26) 
67 ^ 

with the same sign convention, hoc is given by a similarly 
modified equation: 
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With a warped bulk the constraint from hi (7 < 1/12 
in the Minkowski case) is modified. In our equation for 
hi we effectively have two bounds from the two square 
root terms. In the = case only the inner term is of 
any consequence (giving rise to the bound 7 < 1/12). 
When 7^ there are two bounds which are applicable 
in different regimes. If we consider the bound from the 
inner square root we get: 



7< 
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This is valid for > —2. Considering the outer square 
root term we get the bound: 



7< 
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(29) 



This bound becomes negative when > — 1 which is 
unallowed due to our string constraint. The bound on 7 is 
given by the lower of the two constraints when — 2 < < 
— 1 i.e. when in the range where both exists. Therefore 
the bound on 7 for hi to be real is given by (see Fig. llO|l : 




FIG. 10: In order for hi to be real 7 must lie beneath the solid 
line (7m). The region to the right of the vertical dotted line 
is where GBIG4 {he > 0) is allowed (except where = 0). 



we can split the 7, plane into three sections, see Fie. 1101 
The solid line is constructed from the bounds in Eq. H3U|I . 
For hi to be real we must choose values below this line. 
The area to the right of the vertical dotted line (but ex- 
cluding = 0) allows GBIG4. Points to left of this line 
have GBIGl collapsing after a minimum energy density 
(/ib) is reached. The dotted curve on the left comes from 
initial constraint 7 < —4/30. 

In order to consider the a, 7 plane as we did in the = 
Minkowski case, we need the modified Ci^e equations: 



1 - I87 + 27720 ± (1 - 127 - 18720)3/2 

54^ ' 



(32) 



We now have two more formulas for the collapse density 
of GBIG3 and the bounce density of GBIG4: 
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For a particular value of the maximum value of 7 
allowed by this constraint is denoted 7^. When < —9 
the constraint in Eq. H12|l is tighter than that in Eq. I|29() 
(7M < 7m)- This means that for < — 9, /ii is always 
real for allowed values of 7. 

There is a bound for GBIG4 to exist, found by consid- 
ering he = 0. This boimd is given by: 



(31) 



This is a solution to the quadratic obtained from /i^ = 0, 
the other root of the quadratic does not obey the con- 
straint < —4/37 so is ignored. Using our constraints 



(33) 



with the plus sign corresponding to the collapse and the 
minus to the bounce. 

We shall consider the cr, 7 plane for three different val- 
ues of corresponding to three distinct regions in Fig. 1101 
We shall first consider = —1/2. 



1. Typical example = —1/2 

If we take = — 1 /2 , we are in the region where GBIG4 
is allowed. The a, 7 plane can be seen in Fig. ^2 The 
regions I, II and III in Fig.^]extend up to 7m = 8/3. 

In each region we have the following cosmologies. 

• / : (Tc(7m) < cr < 0, 7 > 7i and a < cro(7„i), 7 > 
7m. GBIGl-2 do not exist. GBIG3 reaches a mini- 
mum energy density (/ic) and then evolves back to 
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FIG. 11: The (a, 7) plane for solutions in a AdS {(p = -1/2) 
bulk. The short dotted horizontal line is the maximum value 
of 7 as obtained from Eq. Hl'Zl . 



ji = 00. GBIG4 starts at (0, — /loo) then bounces at 
(/ib, 0) before evolving to (0, +/ioo)- When cr = Cb 
GBIG4 exists as a Minkowski universe. 

• // : cTb < cr < CTc, 7 > 7i- GBIGl-2 and 4 do not 
exist. GBIG3 reaches a minimum energy density 
(/Uc) and then evolves back to /i = cx). When a — (Tc 
GBIG3 ends in a Minkowski universe. 

• /// : cr > (Tc, 7 > 7i- GBIGl-2 and 4 do not exist. 
GBIG3 evolves to a vacuum de Sitter universe. 

• IV : cr < (Tc < 0, 7 < 7m- GBIGl evolves to 
{lic,ho). GBIG2 evolves to ft-oo- GBIG3 reaches a 
minimum energy density (/ic) and then evolves back 
to /i = 00. GBIG4 starts at (/ic,— ^c), bounces at 
(/ib,0) before evolving to (/Ze,+/ic)- When cr = 
CTc GBIGl and 3 evolve to (0, /lo). GBIG4 starts 
at (0, — /lo) and bounces before evolving to (0, /ic). 
When 7 = 7m GBIGl ceases to exist {hi = he). 

• V : CTe < cr < CTb, 7 < 7i. GBIGl-2 both end 
in vacuum de Sitter universes with different values 
of /loo- GBIG3 reaches a minimum energy density 
(/ic) and then evolves back to // = c». GBIG4 starts 
at (0, —/loo) then bounces at (/ib, 0) before evolving 
back to (0, +/100)- When cr = CTb GBIG4 exists as 
a Minkowski universe (0, 0)- 

• VI : CTb < cr < ac, 7 < 7i- GBIGl-2 both end 
in vacuum de Sitter universes with different values 
of /loo- GBIG3 reaches a minimum energy density 
{lie) and then evolves back to /i = 00. GBIG4 
does not exist. When ct = CTc GBIG3 ends in a 
Minkowski universe. When 7 = 7i GBIGl-2 live at 
(0,/ii). 





FIG. 12: The GBIG4 solution in region I for <?!> = -1/2,7 = 
0.1 > 7m and a — —2. The bottom plot is h" vs fi for this 
solution. This shows the differences between this solution 
and the similar case when (j) = —4, where GBIG4 is no longer 
allowed. Arrows denote proper time. 



• VII : CT > CTc, 7 < 7i. GBIGl-3 all end in vac- 
uum de Sitter universes with different values of . 
GBIG4 does not exist. When 7 = 71 GBIGl-2 live 
at (0,/ii). 

In region I there is a GBIG4 solution, which is modified 
in the same way as the GBIG3 solution in the Minkowski 
bulk. As hi^c are no longer valid in region I GBIG4 joins 
onto GBIG2, see Fig.^l The bottom plot in this figure 
shows h" , in order to clearly distinguish this solution 
from the one in the (f) = —4 case where GBIG4 is no 
longer allowed. 

In Fig. ^1 are the results for hoo with </> = —1/2. The 
thin-dark line (ct = —2) lies in regions I and IV. Therefore 
we have only one solution for /loo, which corresponds to 
GBIG2 for 7 < 7ni. For 7m < 7 < 7m this root now 
corresponds to the end point for GBIG3 as in the (j) — 
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FIG. 13: /too for solutions in a AdS {cj> = -1/2) bulk. The 
thin-dark line has a = —2, thin-hght hues have a = —1.5, 
thick-dark line has a — and the thick-light lines have a — 2. 

case. The thin-light hues {a = —1.5) lie in regions I and 

V. There are the two GBIGl-2 solutions and the GBIG4 
solution which converges with the light-thick line at the 
bottom. The thick-dark line {a = 0) lies in regions II and 

VI, so only has GBIGl-2 present. The thick-light lines 
(cr = 2) lie in regions III and VII, so has GBIGl-2 and 
the GBIG3 solution (the horizontal line). 

2. Typical example: (f> — —4 

We now consider ((> — —4 which is in the region 
where GBIG4 no longer exists (see Fig. I14|) . as /ic = 0. 
Therefore fie is not relevant and GBIGl bounces at /ib. 
This means that the (cr, </)) plane. Fig. 1151 is simpler. 
The regions in Fig. 1151 are: 

• I : a < (Tb, 7m < 7 < 7m- GBIGl does not exist. 
GBIG2 starts at (0, — /loo), collapses to (/ib, 0) and 
then expands back to {0,hoo)- GBIG3 expands to 
(/ic,0) and then collapses. When a = GBIG2 
exists as a Minkowski universe. 

• II : (Tb < cr < CTc, 7i < 7 < 7m- GBIGl-2 do not 
exist. GBIG3 expands to /ic and then collapses. 
When a ^ (Jc GBIG3 evolves to a Minkowski uni- 
verse. 

• /// : cr > crc, 7i < 7 < 7m- GBIGl-2 do not exist. 
GBIG3 evolves to ft-oo- 

• IV : cr < crb, 7 < 7ni. GBIGl evolves from /ij to 
/ib and then expands back to /ii. GBIG2 evolves 
to hoo- GBIG3 expands to Hc and then collapses. 
When 7 = 7m GBIGl ceases to exist, GBIG2 ex- 
pands from (/ib,0). When cr = crb GBIGl evolves 
to a Minkowski universe. 
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FIG. 14: Solutions of the Friedmann equation (h vs /x) with 
negative brane tension [a = —4) in an AdS bulk ((^ = —4) 
with 7 = 1/20. For this value of (/> GBIG4 no longer exists 
and GBIGl can collapse. The curves are independent of the 
equation of state w. The arrows indicate the direction of 
proper time on the brane. 



7m 



I 


II 


\ 
\ 

\ 

III 


7m / 




\ 
\ 
\ 


IV /( 






/^b(7) 


V 


\vi 


-4 -2 




a 


2 4 



FIG. 15: The (ct, 7) plane for solutions in a AdS ((^ = —4) 
bulk. The short dotted horizontal line is the maximum value 
of 7 as obtained from Eq. 12811 . The top horizontal line is 
from the initial bound in Eq. UlUII . 



• V : CTb < cr < crc, 7 < 7i. GBIGl-2 evolve to /loo- 
GBIG3 expands to /^c and then collapses. When 
7 = 7i GBIG 1-2 exist as the same dc Sitter uni- 
verse with (0, ft-oo)- When cr = Cc GBIG3 ends as a 
Minkowski universe. 

• VI : cr > crc, 7 < 7i. GBIGl-3 evolve to haa- 
When 7 = 7i GBIGl-2 live at (0, h^). 

In region I we again have a combined solution as 
GBIGl has vanished. As GBIG4 is not allowed (/lo = 0) 
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FIG. 17: hoc for solutions in a AdS (<?!) = -4) bulk. The thin- 
dark line has a = —4.4, thin-light lines have a — —3.2, thick- 
dark line has cr = and the thick-light lines have a — 4.4. 



FIG. 16: The GBIG4 solution in region I for = -4, 7 = 
0.2 > 7m and a — —4. The bottom plot is h" vs fj, for this 
solution. This shows the differences between this solution and 
the case when — —1/2. Arrows denote proper time. 



when we take 7 > 7ni, which causes hi = 0, GBIG2 
matches up with its negative counterpart. We see the 
bouncing GBIG2 solution in Fig. El Note the nature of 
h" is very different to that of the GBIG2, 4 bounce in 
the (j) ~ —1/2 case. In Fig.El^c present results for h^o 
for (h = -4. 
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The (cr, 7) plane for solutions 


in a AdS ((/) = 



bulk. The top horizontal line is the bound from our initial 
bound in Eq. IfT^ . 



to fic and then collapses. When cr = ab GBIGl 
ends in a Minkowski universe. 



3. Typical example: tj) = —20 

Here we shall present results for = — 20 for complete- 
ness. This value of </> lives in the region of Fig. 1 101 where 
hi is always real. This means that the a, 7 plane is much 
simpler. Fig. ITHl 

The regions in Fig. 1181 are: 

• / : cr < CTb, 7 < 7m- GBIGl evolves from fii to 
and back. GBIG2 evolves to ft-oo- GBIG3 expands 



• // : cri < cr < o-c, 7i < 7 < 7m. GBlGl-2 do not 
exist. GBIG3 expands to /ic and then collapses. 
When CT = crc GBIG3 evolves to a Minkowski uni- 
verse. 

• /// : cr > cTc, 7i < 7 < 7m. GBIGl-2 do not exist. 
GBIG3 evolves to hoo- 

• IV : CTb < cr < CTi and crb < cr < ctc; 7 < 7m and 
7 < 7i- GBlGl-2 do not exist. GB1G3 expands to 
fic and then collapses. When 7 = 71 (and for a = 
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FIG. 19: hoc for solutions in a AdS (<^ = -20) bulk. The 
thin-dark line has a = —10, thin- light lines have a = —2, 
thick-dark lines has cr = 2 and the thick-light lines has a = 10. 

(Ti) GBIGl-2 both live at {0,hoo)- When a = CTc 
GBIG3 evolves to a Minkowski universe. 

• V : CT > fTc, 7 < 7i- GBIGl-3 evolve to /loo- When 
7 = 7i GBIGl-2 both live at (0, h^). 

As we decrease the value of 0, region II in Fig. ITHl 
shrinks (the point cri(7M) becomes increasingly positive). 
For (j) < —256/9 region II no longer exists. The nature 
of the solutions in the other regions are unaffected. 

There are no bouncing solutions in this case (in fact 
for any case with cj) < —9) as GBIG4 is unallowed and 
7m < 7m (which rules out the combined solutions). In 
Fig. El we show the hoo results for (j) = —20. 

4. CONCLUSIONS 

In this work we have looked at the general GBIG 
model. We have seen that there is a range of possible 
dynamics that can be achieved depending on the param- 
eters in the model. Our main model of interest is still 
GBIGl as this is the one that starts in a finite den- 
sity "quiescent" singularity and evolves to a de Sitter 
universe. If we warp the bulk enough, GBIGl can be 
allowed to collapse back to its initial density, provided 



there is sufficient negative brane tension. If the bulk is 
warped enough to allow GBIGl to collapse this means 
that the bouncing cosmology GBIG4 can no longer ex- 
ist. The exact nature of the late time dynamics can be 
changed by including a non-zero brane tension. A neg- 
ative brane tension can reduce the Hubble rate at late 
time and a positive tension will increase it. GBIGl can 
end in a future "quiescent" singularity with a non-zero 
density if we have an appropriate (negative) brane ten- 
sion present. As the solution of interest is GBIGl and 
we want it to provide the late time acceleration that we 
are experiencing it needs to end as a vacuum de Sitter 
universe. Therefore if there is some brane tension it must 
take values Ci > cr > CTc- 

GBIG2 also starts in a "quiescent" singularity but it 
super-accelerates so is un-physical and of little interest. 

GBIG3 still starts in an infinite density big-bang but 
has a number of possible late time dynamics. In a 
Minkowski bulk GBIG3 can either evolve to a Minkowski 
state, as shown in Ref. a vacuum de Sitter state or 
even loiter around fi = —a before ending in a vacuum de 
Sitter state or a "quiescent" singularity. This loitering 
cosmology is different from that in Ref. ^^I; ^ do 
not require a naked bulk singularity or a de Sitter bulk. 
If we warp the bulk GBIG3, will generally collapse, un- 
less there is sufHcient (positive) brane tension to allow 
the solution to end in a vacuum de Sitter state. 

GBIG4 can only ever exist in a mildly warped bulk 
with negative brane tension. So it can be said that 
GBIG4 is un-physical due to requirement that ct < 0. 

There are a number of bouncing cosmologies within 
this set-up. There are the GBIG4 solutions as mentioned 
above. There are also the solutions where GBIG4 and 2 
match up (7m < 7 < 7m with (/)GBiG4Lim < (f> < 0) and 
the GBIG2 bouncing cosmologies (7,„ < 7 < 7m with 
—9<4>< 0GBiG4Lim)- Each of these have different dy- 
namics so would produce different evolutionary histories. 
The solutions that spend time on the GBIG2 branch will 
experience phantom like behaviour during this period. 
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